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ABSTRACT 

Laminar  /  turbulent  transition  is  studied  in  a  straight  channel  with  2  Hz  imposed 
bulk  flow  unsteadiness  for  Reynolds  numbers  of  1 100  to  3800  and  Strouhal  numbers  of 
0.121  to  0.035.  Channel  aspect  ratio  is  40  to  1  with  1.27  cm  height,  50.8  cm  width  and 
4.27  m  length.  Observations,  videos,  and  photographs  of  smoke  patterns  show  different 
subcritical  transition  events  including  Tollmien-Schlichting  waves,  vortex-array  type 
motion  evidenced  by  smoke  swirls,  ribbon-like  patterns,  turbulentspots,  and  fully  tur- 
bulent flow.  Fluctuating  intensity  magnitudes  relative  to  velocities  from  imposed  un- 
steadiness are  determined  from  phase-averaged  velocity  traces.  With  2  Hz  unsteadiness, 
transition  events  begin  to  occur  at  lower  Reynolds  numbers  and  extend  over  a  wider 
range  of  Reynolds  numbers  compared  to  flows  with  no  imposed  unsteadiness. 
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L   INTRODUCTION 

A.     BACKGROUND 

Many  typical  flows  found  in  nature  and  in  aerodynamic  applications  and  almost 
all  flows  associated  with  turbomachinery  involve  freestream  unsteadiness.    Most  of 
these  flows,  especially  those  of  technological  interest  undergo  transition  from  a 
laminar  state  to  a  turbulent  condition.  The  accurate  prediction  of  the  location  and 
extent  of  transition  and  the  accompanying  changes  of  important  flow  properties  is 
vital  to  the  design  of  aircraft  and  most  energy  transformation  machines  involving 
fluid  flow. 

Laminar-turbulent  transition  results  in  significant  increases  in  local  wall  shear 
stress  and  convective  heat  transfer  rates.  This  becomes  evident  by  first  considering 
laminar  flows,  where  transport  of  heat  and  momentum  normal  to  streamlines  occurs 
only  as  a  result  of  molecular  diffusion.  This  is  usually  an  advantage  since  tangential 
forces  at  solid  surfaces  (drag)  and  heat  flux  normal  to  surfaces  are  relatively  low. 
Occasionally,  laminar  flows  are  undesirable  since  they  tend  to  separate  from  surfaces 
more  easily  than  turbulent  flows.   In  most  applications,  laminar  flows  are  unstable  and 
transitional  or  fully  turbulent  flows  exist.  In  the  latter,  higher  rates  of  diffusion  across 
mean  streamlines  due  to  turbulent  transport  cause  significantly  higher  drag  and  heat 
flux  at  surfaces. 

The  parameter  generally  used  to  describe  the  onset  of  transition  is  the  Reynolds 
number.    In  channel  flows,  linear  stability  theory  predicts  that  the  laminar  flow 
becomes  unstable  to  small  amplitude  disturbances  at  a  critical  Reynolds  number  of 
5772.    Here  the  Reynolds  number  is  based  on  mean  centerline  velocity  and  channel 
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half- width.   With  large  amplitude  disturbances,  however,  the  transitional  Reynolds 
number  for  channel  flow  may  be  as  low  as  1000-1500.  The  value  of  Reynolds  number  at 
the  onset  of  transition  depends  strongly  on  the  disturbances  existing  in  the  flow  and 
conditions  which  prevail  in  the  entry  region  of  the  channel. 

Also,  the  sequence  of  events  which  occur  during  transition  are  strongly  dependent 
on  upstream  history,  especially  laminar  boundary  layer  character  and  development. 
There  is  a  growing  collection  of  evidence  which  also  suggests  that  turbulent  shear  layer 
development  is  more  strongly  dependent  on  its  upstream  history  than  previously 
believed.   This  would  mean:  (1)  that  transition  affects  turbulent  boundary  layer 
structure  downstream,  and  (2)  that  turbulent  boundary  layer  results  obtained  under 
equilibrium  conditions  (i.e.  streamwise  self-similarity)  may  not  adequately  represent 
flow  behavior  near  surfaces  of  operating  turbines  and  airfoils.    Non- equilibrium  shear 
layer  development  near  turbine  passage  surfaces  and  near  airfoils  is  due  to  varying 
surface  curvature,  varying  pressure  gradients,  and  freestream  unsteadiness. 
Unsteadiness  is  particularly  important  in  this  regard  since  it  affects  shear  properties  at 
a  given  location  in  two  ways:  (1)  by  its  local  direct  influence  and  (2)  from  its  affect  on 
shear  layers  upstream,  whose  history  affects  downstream  development.    The 
computational  work  of  Singer,  Ferziger  and  Reed  [Ref.  1]  indicates  that  moderately  low 
frequency  oscillations  in  the  channel  can  initiate  nonlinear  effects  which  trigger 
transition  on  flow  regimes  where  the  steady  flow  remains  laminar.   These  happen  even 
though  linear  theory  indicates  that  the  flow  is  stable.    Experimental  clarification  of 
these  points  is  needed,  along  with  the  verification  of  the  sequence  of  events  occurring 
during  the  transition  process  at  different  Strouhal  numbers.  The  view  of  the  overall 
transition  process  is  especially  valuable  to  the  gas  turbine  research  community. 


Unsteadiness  is  additionally  important  since  boundary  layers  which  exist  near 
turbulent  blades  are  significantly  different  from  equilibrium  boundary  studied  in 
many  laboratories.  A  number  of  recent  studies  indicate  the  inadequacy  of  results 
obtained  from  equilibrium  boundary  layers.   Gillis  and  Johnston's  [Ref.  2] 
measurements  in  boundary  layers  downstream  of  a  convex  curved  surface  and  Ligrani's 
[Ref.  3]  results  in  boundary  layers  downstream  of  artificial  thickening  devices  showed 
that  many  structural  changes  persist  long  distances  downstream  of  perturbations.  This 
indicates  that  boundary  layer  behavior  may  depend  on  upstream  history  to  a  greater 
extent  than  previously  believed.   This  is  important  regarding  unsteadiness  since  it 
alters  boundary  layers  locally  as  well  as  the  history  of  boundary  layers  upstream. 

Few  experiments  have  been  performed  to  study  the  stability  or  transition  of 
oscillating  flow  in  a  plane  channel.    Most  known  information  is  from  a  handful  of 
analytic  studies  and  the  full  Navier-Stokes  simulation  undertaken  at  Stanford. 

Regarding  the  latter,  the  work  by  Singer,  Reed  and  Ferziger  [Ref.  1]  provides 
fascinating  information  as  to  the  effects  of  unsteadiness  on  transition.   For  plane 
channel  flow  linear  theory  has  been  verified  at  a  Reynolds  number  of  5000.  With  mean 
flow  frequency  in  the  range  of  the  non-dimensional  frequency  of  the  least  stable  TS 
wave  (Tollmien-Schlichting),  they  find  that  the  mechanisms  driving  the  transition  are 
similar  to  those  described  by  Herbert's  secondary  instability  theory  [Ref.  4  and  51  and 
observed  experimentally  by  Nishioka  et.  al.  [Ref.  6  and  7],  and  Kozlow  and  Ramazanov 
[Ref.  81.  In  the  plane  channel  flow  these  mechanisms  have  been  studied  extensively  in 
previous  computational  investigations.   The  addition  of  mean  flow  modulation  causes 
significant  variations  in  the  amplitude  of  the  TS  wave  during  the  period,  hence  the 
dominant  flow  structures  depend  on  the  Strouhal  number  and  the  magnitude  of  the 
modulation. 


In  a  straight  channel  without  imposed  unsteadiness  the  process  of  transition 
begins  with  TS  waves  experiencing  mild  growth  at  sufficiently  high  Reynolds  numbers. 
When  these  waves  grow  to  an  amplitude  of  about  0.3  percent  of  the  free  stream  velocity, 
three-dimensional  waves  become  unstable.  Transition  often  occurs  at  Reynolds 
numbers  below  that  at  which  TS  waves  become  unstable.  These  waves  cause  Lambda 
shaped  vortices  to  lift  from  the  near  wall  region  and  "burst"  ejecting  low  speed  fluid 
from  close  to  the  wall  into  the  faster  moving  fluid  away  from  the  wall.  These  events 
then  lead  to  a  rapid  breakdown  of  the  laminar  boundary  layer  into  a  turbulent  one.  The 
process  can  be  modified  by  a  number  of  different  factors.   Freestream  turbulence  can 
lead  to  non-sinusoidal  oscillations  in  the  nominally  laminar  boundary  layer.    Local 
disturbances,  such  as  roughness  elements,  may  lead  to  structures  known  as  turbulent 
spots.  These  patches  of  turbulent-like  flow  form  a  characteristic  arrowhead  shape  and 
spread  in  the  downstream  direction.  Distributed  roughness,  acoustic  disturbances, 
surface  waviness  and  other  factors  all  serve  to  complicate  the  basic  process.  There  also 
exist  equilibrium   states  where  two  dimensional  waves  (which  are  similar  in  form  to  TS 
waves)  with  a  given  amplitude  will  neither  grow  nor  decay.   The  initial  disturbances  in 
the  inlet  of  the  channel  can  evolve  into  these  or  similar  structures,  after  which  the 
strong  secondary  instability  of  three-dimensional  waves  can  take  place. 

B.      OBJECTIVES 

This  research  is  a  follow-on  to  the  work  undertaken  by  Longest  |Ref.  9],  and  is 
believed  to  be  the  first  experimental  (quantitative)  study,  known  to  the  author,  on  the 
effects  of  imposed  mean  flow  unsteadiness  on  transition  in  straight  channel  flow.  The 
objective  is  to  investigate  the  effects  of  imposed,  periodic,  and  deterministic 


unsteadiness  on  the  transition  process  from  laminar  to  turbulent  flow.  Answers  to  the 
following  questions  are  sought: 

1.)  What  parameters  govern  the  unsteady  flow? 

2.)  What  events  occur  in  the  flow  during  the  onset  and  development  of 
transition? 

3.)  What  is  the  interaction  between  the  organized  induced  unsteadiness  and 
transition  (e.g.,  development  of  two-dimensional  Tollmien-Schlichting  waves  into 
three-dimensional  instabilities)? 

4.)  At  what  range  of  frequencies  can  one  expect  interactions  between  the 
organized  imposed  unsteadiness  and  transition  phenomena? 

Flow  measurements  (using  hot-wire  anemometry)  and  flow  visualization  are  used 
to  analyze  the  effect  of  oscillating  the  bulk  flow  in  otherwise  steady  plane  Poiseuille 
flow.  An  attempt  is  made  to  determine  whether  the  type  of  unsteadiness  studied 
(sinusoidal  at    two  hertz)  contributes  to  the  early  onset  of  transition,  delay  of 
transition,  or  has  no  effect.   Sinusoidal  oscillation  was  chosen  to  induce  a 
deterministic  periodic  variation  of  velocity.  The  results  are  compared  with  the  results 
of  Singer's  [Ref.  1]  numerical  simulation,  which  suggests  that  such  oscillations  provide 
a  stabilizing  effect  at  all  but  very  low  frequencies. 

C.     ORGANIZATION 

Section  II  follows  this  introduction  and  describes  the  experimental  facilities.    In 
section  III  experimental  procedures  are  discussed,  including  calibrations,  channel 
validation,  flow  measurements  and  flow  visualization.   The  phase  averaging  technique 
is  also  detailed  here.  A  discussion  of  the  results  follows  in  section  IV.   Conclusions  are 
presented  in  section  V. 


H.   EXPERIMENTAL  FACILITIES 

A.       CHANNEL 

The  40  to  1  aspect  ratio  straight  channel,  shown  schematically  in  Figure  1,  is 
located  in  the  laboratories  of  the  Department  of  Mechanical  Engineering.  The  details 
of  the  design  and  construction  are  discussed  by  Ligrani  and  Subramanian  [Ref  10]. 
Briefly,  the  channel  is  made  of  6.35  mm  thick  plexiglass  with  a  straight  test  section 
4.27  m  in  length,  with  inside  dimensions  of  1.27  cm  in  height  and  50.8  cm  in  width.  The 
top  and  bottom  walls  of  the  channel  are  supported  by  ribs  and  cross  beams  along  its 
length.  The  side  walls  are  removable  in  order  to  gain  access  to  the  inside  of  the  channel. 
The  inlet  section  consists  of  a  honeycomb  assembly,  3  frames  with  screens  and  a  20: 1 
contraction  ratio  nozzle.   Two  layers  of  cheesecloth  are  placed  at  the  inlet  lip  for 
filtering  the  air  entering  the  channel.  At  the  exit  of  the  4.27  m  long  test  section  are 
three  10. 16  cm  long  frames  with  screens  for  housing  the  unsteady  device  (the  device  used 
to  create  unsteadiness  in  the  flow)  and  a  second  honeycomb.  A  two  dimensional 
diffuser  45.72  cm  long  with  a  3  degree  total  angle  is  located  just  upstream  of  the  plenum 
chamber.  The  exit  plenum  chamber  has  inside  dimensions  of  60.96  cm  by  60.96  cm.  A 
3.81  cm  orifice  plate  assembly  between  the  plenum  chamber  and  the  blower  is  provided 
to  meter  the  flow  through  the  channel.  The  exit  of  the  plenum  chamber  is  connected  to 
the  suction  side  of  a  5  H.P.  blower  with  a  long  50.8  cm  diameter  pipe.  A  valve  is  provided 
between  the  plenum  and  the  blower  to  regulate  the  flow.  Bulk  flow  velocities  up  to  15 
m/s  are  achieved  in  the  channel.    In  our  experiments,  bulk  flow  velocities  in  the  range 
1.0-4.5  m/s  are  used. 


The  channel  is  designed  with  a  sufficient  length  such  that  transition  occurs  after 
the  laminar  velocity  profile  (parabolic)  is  fully  developed  (self  similar).    The  initial 
flow  development  length  is  estimated  (Ref.  1  l:p.  186)  to  be  within  about  25  per  cent  of 
the  test  section  length  for  the  range  of  speeds  used  in  the  present  tests. 

B.     UNSTEADY  DEVICE 

The  most  important  aspect  in  the  design  of  any  unsteady  device  for  a  transition 
study  is  to  ensure  that  the  device  produces  a  deterministic  and  periodic  unsteadiness 
without  adding  other  disturbances.  To  achieve  this  in  an  open  circuit  induction  tunnel, 
the  unsteady  device  is  best  located  at  the  exit  of  the  channel.  This  way,  problems  caused 
from  the  wakes  and  the  flow  disturbances  resulting  when  the  unsteady  device  is 
upstream  of  the  test  section  are  avoided. 

The  design  is  based  on  the  design  used  by  Miller  and  Fejer  [Ref.  12]  and  is  shown  in 
Figure  2.   Since  the  depth  of  the  channel  is  only  12.7  mm,  a  single  rotating  vane  driven 
through  a  spur  gear  train  by  a  DC  stepper  motor  is  adequate  to  introduce  the  required 
unsteadiness  in  the  flow.  The  rotating  vane  causes  varying  resistances  to  the  flow  in 
the  test  section.  The  frequency  of  the  imposed  unsteadiness  is  controlled  by  changing 
the  vane  rpm  and  the  amplitude  is  controlled  by  using  different  vane  widths.   It  is  also 
important  to  match  the  flow  speed  with  the  vane  speed.  A  vane  8.7  mm  wide  was  used  in 
this  study. 

The  vane  in  Figure  2  is  made  of  a  3.2  mm  thick  brass  strip  with  rounded  edges,  and 
spans  the  entire  width  of  the  channel.  The  vane  is  supported  at  the  ends  by  a  3.2  mm 
diameter  shaft  and  bushings  that  are  fitted  to  the  side  walls  of  the  frame.  Three 
additional  intermediate  spanwise  struts  are  provided  to  increase  the  rigidity  of  the 
vane.  One  end  of  the  shaft  is  extended  to  accommodate  a  48  TPI  (Threads  Per  Inch)  spur 


gear.  The  driver  spur  gear  has  12  TPI  and  is  mounted  on  the  Superior  Electric,  M092- 
FD310  Stepper  Motor.  To  achieve  a  fine  step  control,  a  Modulynx  MITAS  PMS085- 
D050  Drive  with  Modulynx  MITAS  PMS085-C2AR  Drive  Controller  is  used.   This 
combination  facilitates  motor  shaft  speeds  up  to  100  revolutions  per  second,  with  a 
constant  200  steps  per  motor  shaft  revolution.  Note  that,  with  a  1:4  gearing  ratio  and 
two  cycles  of  imposed  flow  oscillations  per  vane  rotation,  one  motor  revolution 
corresponds  to  one  half  cycle  of  the  imposed  flow  oscillation  by  the  vane.  The  Drive 
Controller  is  programmed  for  the  required  motor  speed  (rps).  The  Modulynx  MITAS 
PMS085-D050  Drive  was  also  used  to  provide  a  trigger  signal  to  initiate  data 
acquisition.   This  trigger  also  enabled  the  correlation  of  data  taken  with  the  position  of 
the  vane  at  the  time  that  the  data  was  recorded. 

C.       FLOW  MEASUREMENTS 

The  equipment  used  for  flow  measurement  includes  a  digital  manometer,  hot-wire 
probe,   hot-wire  bridge,  signal  conditioner,  digital  multimeter,  high  speed  data 
acquisition  system,  oscilloscope  and  a  microcomputer.   These  are  described  in  detail 
below 

1.  Bulk  Velocity 

The  bulk  flow  velocity  through  the  channel  is  determined  from  the  pressure  drop 
across  the  3.81  cm  orifice  plate  which  is  installed  in  the  piping  between  the  blower  and 
the  exit  plenum  chamber.  This  pressure  drop  is  measured  and  displayed  by  a  Validyne 
Model  PS  309  digital  manometer  (inches  of  water). 


8 


2.    Instantaneous  Velocity 

A  constant  temperature  hot-wire  anemometer  is  used  to  measure  instantaneous 
velocities  of  the  flow.    Briefly,  the  principles  of  hot-wire  anemometry  (constant 
temperature  mode)  are  as  follows. 

The  flow  sensing  element  (hot-wire  probe)  consists  of  a  short  length  of  wire, 
connected  to  two  supports,  which  forms  the  leg  of  a  wheatstone  bridge.  This  bridge 
maintains  the  wire  at  elevated  temperature  and  is  initially  balanced  prior  to  exposing 
the  probe  to  the  flow.  This  balancing  resistance,  referred  to  as  "cold"  resistance  is 
multiplied  by  an  "overheat  ratio".  Once  the  bridge  is  set  to  this  value,  and  as  the  probe  is 
exposed  to  the  flow,  the  bridge  is  placed  in  operation  and  a  small  current  (60  mA)  passes 
through  the  circuit.  Wire  temperature  increases  because  of  resistive  heating  and  the 
wire  resistance  increases  until  the  bridge  is  rebalanced.  The  governing  relationship  is: 

R=Ro(l  +a(T-To)) 
where  R  is  the  resistance  at  temperature  T,  Ro  is  the  resistance  at  the  reference 
temperature  To,  and  ais  the  temperature  coefficient  of  resistance. 

Flow  past  the  probe  causes  the  wire  to  cool  (by  convection)  and  the  bridge  becomes 
unbalanced  as  the  resistance  decreases.  The  bridge  is  rebalanced  by  a  feedback  loop  and 
amplifier.   Response  is  rapid  enough  to  maintain  wire  temperature  and  resistance 
nearly  constant.   As  such,  the  variations  in  the  bridge  voltage  are  directly  proportional 
to  the  current  variations,  and  flow  velocity  can  be  calculated  as  a  function  of  the  bridge 
voltage. 

Now  the  different  components  used  to  measure  instantaneous  velocity  are  described. 

a.  Hot-wire  Probe  and  Probe  Positioning 

A  Dantec  Electronics  Inc.  P51  probe,  mounted  such  that  the  wire  is  horizontal  and 
normal  to  the  streamwise  flow  direction,  is  used.  The  wire  sensor  is  5.0  microns  in 


diameter  and  approximately  1.25  mm  long.  The  probe  is  mounted  through  the  side  wall 
of  the  channel  and  is  held  fixed  by  a  rotatable  lever  arm  which  allows  placing  of  the 
probe  at  various  y/d  positions  ranging  from  0.2  to  0.9.  A  scale  graduated  in  terms  of  y/d 
is  mounted  adjacent  to  the  lever  arm  for  probe  positioning  to  an  accuracy  of 
approximately  0.5  mm.   Longitudinal  and  transverse  position  of  the  probe  for  this 
study  were  46. 12  cm  and  15.64  cm,  measured  from  the  near  flange  at  the  end  of  the  test 
section,  and  from  the  inside  of  the  left  channel  wall  (looking  downstream), 
respectively.  By  rearranging  the  sidewalls  the  probe  can  be  positioned  at  other 
streamwise  locations  for  future  studies. 
b.  Hot-wire  Bridge 

The  hot-wire  is  operated  with  a  DISA  55  M10  constant  temperature  bridge.  A  five 
meter  coaxial  cable  is  used  to  connect  the  hot-wire  to  the  bridge.  For  the  purposes  of 
this  study  the  cold  resistance  was  multiplied  by  an  overheat  ratio  of  1.8.  The  DC  voltage 
from  the  bridge  is  measured  with  a  Hewlett-Packard  3466A  digital  multimeter. 
c  Signal  Conditioner 

The  hot-wire  signal  is  filtered  and  amplified  by  a  DISA  56W20  signal  conditioner. 
The  high-pass  filter  is  set  to  0. 1  Hz  in  all  cases  to  remove  the  DC  signal.  This  however, 
did  not  alter  the  effects  of  2.0  Hz  imposed  flow  unsteadiness.    The  low-pass  filter  is  set 
at  10  KHz  to  filter  out  noise  existing  above  this  frequency.  The  hot-wire  signal  is 
amplified  using  a  gain  of  10  to  achieve  an  acceptable  signal  to  noise  ratio.   Output  from 
the  signal  conditioner  is  sent  to  the  high  speed  data  acquisition  system  and  to  a  B  &  K 
Precision  10  MHz  Oscilloscope  used  to  signal  monitor  the  signal. 
d.  Highspeed  Data  Acquisition  System 

A  Hewlett-Packard  6944A  Multiprogrammer  is  used  in  conjunction  with  a 
Hewlett-Packard  series  9000,  Model  310  microcomputer  to  acquire  data.  A  buffered 
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69759A  A/D  card  configuration  is  used  with  the  Multiprogrammer  to  allow  continuous 
data  acquisition.  The  A/D  card  is  externally  triggered  using  2-3  volt  TTL  type  signals 
from  the  unsteady  device  motor  drive.  When  the  unsteady  device  is  not  in  operation, 
triggering  is  provided  by  a  Hewlett-Packard  33 1 1A  Function  Generator  set  to  the  same 
frequency  as  that  used  with  the  unsteady  device  in  operation.  Sampling  occurs  as 
square  wave  voltage  changes  from  high  to  low. 

The  storage  capacity  of  the  A/D  buffer  is  64.000  bytes.  The  number  of  imposed 
unsteadiness  cycles  which  can  be  sampled  depends  on  the  desired  number  of  data  points 
per  cycle.  For  example,  if  100  points  per  cycle  are  desired,  then  640  cycles  can  be 
sampled.  In  this  study,  sampling  was  done  at  100  points  per  cycle  for  all  cases,  where 
cycle  refers  to  the  cycle  of  flow  blockage  caused  by  the  vane  rotation  of  the  unsteady 
device.     100  points  per  cycle  were  chosen  so  that  there  would  be  sufficient  resolution  to 
determine  the  effect  of  unsteadiness  over  the  complete  phase  of  each  unsteadiness  cycle. 
e.  Datastorage 

The  Hewlett-Packard  microcomputer  is  used  to  acquire  the  data  from  the  A/D 
buffer,  using  a  data  acquisition  program  HOTWIREPAV,  described  in  the  software 
directory  (APPENDIX  A).   Raw  instantaneous  voltage  data  is  acquired  and  stored  on 
standard  3.5  inch  micro  diskettes  for  further  processing.  Manual  inputs  to  the  program 
include  bulk  flow  velocity,  flow  blockage  frequency,  and  date  and  time. 
f.  Data  Processing 

Processing  of  the  data  is  done  using  a  Hewlett-Packard  microcomputer  using  the 
programs  MAX_MIN  and  AVGPLOT.   This  software  is  described  in  detail  in  the  software 
directory  (APPENDIX  A).   Many  of  the  curves  used  to  graphically  display  the  results 
were  generated  using  a  Macintosh  SE  personal  computer  and  the  graphics  program 
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CRICKETGRAPH.  The  phase  averaging  technique  employed  by  the  program  MAX_MIN 
is  described  below. 

In  the  case  of  an  unsteady  flow  the  instantaneous  velocity  can  be  split  into  three 
components. 

U  =  U  +    U  +    U' 
where  U   is  the  time  averaged  velocity,  0  is  the  imposed  phase  averaged  velocity  and 
Uf   is  the  instantaneous  fluctuating  component.  Q  is  the  component  due  to  the 
unsteadiness,  and  is  zero  for  steady  bulk  flow.  In  the  present  study  U   and  U   are 
combined  to  represent  the  phase  averaged  velocity  (U  ).  Thus, 

u  =  a  +  u' . 

0    is  determined  by  phase  averaging  U  over,  say  640  cycles  (ncycl),  at  different 
positions  of  the  phase  of  each  oscillation. 
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When  u'  is  a  large  percentage  of  U  then  a  larger  number  of  cycles,  sometimes  up  to 
4640  cycles,  are  used  to  determine  U  .  For  this,  multi-run  data  acquisition  is  made,  with 
the  vane  at  the  same  starting  position  for  each  run.   Only  the  last  580  data  points  (due  to 
sampling  procedure)  of  each  run  are  used  in  the  cumulative  phase  averaging.  Once  U  is 
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determined,  U'   is  calculated  from  il'  =   U  -   b  ,  and  {u')    =  (U  -  0)     .  Therefore, 
the  phase  averaged  root  mean  squared  velocity  is. 
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Time  averaged  velocity,  U  t  and  rms  velocity  V  {u')    are  calculated  by  averaging 
over  all  the  points  (100  x  ncycl)  per  run  and  in  the  case  of  multi-runs,  100  x  ncycl  x  no. 
of  multi-runs. 

ncycl  1  oo 


ncycl  t-'A  00^    v     '    ' 

■f  m  =1  n  =1 


ncycl  100 


^=j74^iM(U(m'n)-D) 

V  J  m  =1  n  =1 


D.      FLOW  VISUALIZATION 

The  equipment  utilized  for  flow  visualization  studies  includes  a  smoke  wire  for 
generating  smoke  in  the  channel  and  video  equipment  to  photograph  the  flow  in  the 
channel.  These  are  described  in  detail  below. 
1.  Smoke  Wire  System 

The  smoke  wire  system  is  a  modified  Vertical  Smoke  Wire  Instrument  made  by 
Flow  Visualization  Systems  of  Bolingbrook,  IL.   The  system  in  its  present  configuration 
consists  of  a  0. 127  mm  diameter  nichrome  wire  that  is  stretched  across  the  width  of  the 
channel  through  access  ports  in  the  channel  side  walls.  The  wire  is  passed  through 
these  access  ports  through  a  thin  stainless  steel  tube  (hypodermic  needle)  which  allows 
the  wire  to  move  freely.  The  smoke  wire  is  connected  to  an  A/C-D/C  converter  and  this 
in  turn  has  voltage  supplied  by  a  Calrad  45-740  (0-130  VAC)  variac  set  at  45  VAC.  The 
wire  is  energized  by  pulsing  the  "Momentary  On"  switch  of  the  A/C-D/C  converter. 
Since  the  wire  acts  as  a  resistor  and  heats  up  when  a  current  is  passed  it  becomes  more 
elastic  and  has  a  tendency  to  sag  in  the  channel.  To  prevent  this  from  affecting  the  flow 
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visualization,  a  tensioning  device  has  been  installed  on  one  end  of  the  wire  to  ensure 
that  the  wire  is  held  taut  in  the  channel  when  it  is  energized.  This  tension  device 
consists  of  two  pulleys  which  support  the  wire  and  a  third  pulley  with  a  weight  attached 
between  the  two  support  pulleys.  The  wire  is  coated  with  a  paraffin  based  oil  which  is 
applied  by  means  of  a  thin  rod  with  an  absorbent  piece  of  felt  at  the  tip.  When  applied  in 
this  manner,  the  wire  is  coated  evenly  with  small  droplets  of  oil.   Energizing  the  wire 
produces  a  thin  even  sheet  of  smoke  until  all  the  oil  on  the  wire  burns.  The  smoke  from 
the  wire  lasts  for  about  4-8  seconds  depending  on  the  flow  rate  in  the  channel. 

The  wire  is  placed  horizontally  normal  to  the  flow  76.8  cm  upstream  from  the 
end  of  the  test  section.   For  the  purposes  of  this  study,  the  wire  was  located  at  a  y/d  value 
of  0.84  (d  being  the  width  of  the  channel  and  y  the  height  from  the  bottom  of  the 
channel). 

2.  Video  Equipment 

The  video  equipment  used  for  flow  visualization  is  the  same  as  that  described  by 
Longest  [Ref.  9].   Briefly,  a  high  speed  video  camera  is  mounted  on  a  tripod  above  the 
channel  and  focused  on  the  area  of  interest.  The  focal  plane  from  the  camera  lens  to  the 
channel  is  surrounded  by  an  enclosure  to  exclude  all  extraneous  light  and  thus 
decreasing  the  reflections  which  interfere  with  the  flow  visualization.   The  bottom  of 
the  channel  was  blacked  out  and  lighting  was  provided  by  two  600  Watt  floodlights 
placed  1  m  from  the  sidewall  on  either  side  of  the  channel  at  the  same  height  as  the 
channel. 

Still  photographs  were  made  by  taking  pictures  of  the  aforementioned  video 
when  viewed  on  a  monitor.  A  Nikon  F-3  SLR  and  Kodak:Tri-Xpan  film  were  used  to 
take  the  pictures.  The  shutter  speed  was  one-eighth  of  a  second  and  the  focal  aperture 
was  set  at  f  5.6. 
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m.  EXPERIMENTAL  PROCEDURES 

A.  ORIFICE  DP  CALIBRATION 

The  orifice  differential  pressure  calibration  for  bulk  velocity  flow  measurements  was 
performed  by  Longest  [Ref.  9).  Results  of  the  calibration  are  shown  in  Figure  4. 

B.  HOT-WIRE  CALIBRATION 

The  hot-wire  was  calibrated  in  the  freestream  of  the  wind  tunnel,  also  located  in  the 
laboratories  of  the  Mechanical  Engineering  Department  of  the  Naval  Postgraduate 
School. 

The  hot-wire  probe  is  mounted  normal  to  the  flow  in  the  center  of  the  channel  so  that 
it  is  in  the  freestream.  The  hot-wire  is  connected  to  the  same  hot-wire  bridge  previously 
discussed.   Output  from  the  hot-wire  bridge  is  read  on  a  digital  multimeter.   Freestream 
velocity  in  the  wind  tunnel  is  measured  utilizing  a  Kiel  pressure  probe,  a  wall  static 
pressure  tap,  and  a  Validyne  digital  manometer.   Voltage  and  differential  pressure 
readings  are  taken  for  a  range  of  pressure  drops  corresponding  to  velocities  of  1.0  to  4.0 
m/s.  The  computer  program  HWCAL  is  run  on  the  Naval  Postgraduate  School's  IBM 
main  frame  computer  to  calculate  all  calibration  constants.    Following  the  calibration, 
the  hot-wire  is  installed  in  the  straight  channel. 

C.  CHANNEL  VALIDATION 
1.  Channel  Flow  Tests 

Assembly  and  installation  of  the  straight  channel  in  the  Mechanical  Engineering 
Laboratory  of  the  Naval  Postgraduate  School,  along  with  the  initial  flow  visualization 
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and  measurement  results  are  detailed  by  Longest  (Ref.  9].  Prior  to  commencing  flow 
measurements  for  this  study  a  check  of  the  channel  flow  quality  is  performed  as 
follows.   Smoke  is  generated  using  the  smoke  generator  described  by  Longest  [Ref.  9]  and 
then  drawn  into  the  channel  at  a  speed  of  approximately  2.0  m/s.  The  first  few  times 
this  was  performed,  numerous  leaks  from  the  sidewalls  of  the  channel  were  observed. 
These  leaks  are  eliminated  using  removable  spring  clamps  on  both  sides  of  the  channel, 
along  its  entire  length.  The  clamps  are  left  in  position  while  the  channel  is  in 
operation. 

2.  Laminar  Average  Velocity  Profile  Verification 

Velocity  profile  verification  is  performed  at  a  Reynolds  number  of  1 103  (bulk 
velocity  of  1.32  m/s)  when  the  flow  is  laminar.   Instantaneous  velocity  measurements 
were  taken  at  eight  y/d  positions  ranging  from  0.2  to  0.9,  at  an  x/d  of  300.0.  This  data 
was  taken  using  the  HOTWIREPAV  program,  with  no  imposed  unsteadiness.   MAX_MIN 
was  used  to  process  the  data,  and  the  time-averaged  local  velocities  were  normalized 
with  respect  to  the  maximum  value  which  occurred  at  the  mid-channel  location  (y/d  = 
0.5).  A  second  order  polynomial  curve  was  fit  to  the  data  and  is  shown  in  Figure  5.  The 
curve  demonstrates  that  the  flow  is  fully  developed  at  this  location  and  follows  the 
expected  parabolic  behavior. 

D.     FLOW  MEASUREMENT 
1.  Reynolds  Number  Survey 

Velocity  surveys  are  conducted  for  Reynolds  numbers  ranging  from  1 100  to  3800 
with  the  hot-wire  probe  at  a  y/d  position  of  0.88,  and  a  z/d  of  -7.7.  This  survey  is  done  to 
determine  how  the  imposed  unsteadiness  affects  the  transition  process  as  Reynolds 
number  varies.   Data  at  each  Reynolds  number  are  taken  by  recording  results  with  no 
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imposed  unsteadiness  followed  by  results  with  imposed  unsteadiness.   Data  sets  of 
64,000  points  are  taken  using  only  the  last  58,000  for  data  processing,  as  previously 
described.  Average  velocity  and  RMS  velocity  values  from  such  survey  are  used  to 
generate  Figures  6  and  7.    These  experiments  were  repeated  at  several  Reynolds 
numbers  to  ensure  consistent  and  repeatable  data  trends. 

2.  Multi-run  Averaging 

Multi-run  phase  averaging  using  up  to  4640  cycles  of  flow  unsteadiness  were 
conducted  at  several  Reynolds  numbers.  The  purpose  of  this  multi-run  averaging  is  to 
show  that  the  survey  data  is  statistically  stable  for  higher  Reynolds  numbers  where 
fluctuation  levels  are  significant  with  respect  to  the  local  mean  velocity.   Results  are 
presented  in  Figures  8-35. 

3.  Velocity  Profiles 

Velocity  profiles  with  and  without  imposed  unsteadiness  are  taken  at  the  same 
Reynolds  numbers  as  used  for  the  multi-run  averaging.  The  purpose  of  these  profiles  is 
to  record  flow  development  as  the  Reynolds  number  is  changed  and  to  study  the  effect  of 
imposed  unsteadiness  on  average  and  RMS  velocities  as  a  function  of  the  Reynolds 
number  and  y/d.   Data  acquisition  and  processing  is  performed  as  described  in  the 
Reynolds  number  survey.   Results  are  provided  in  Figures  36-82. 

E.      FLOW  VISUALIZATION 

Smoke  patterns  were  recorded  on  film  using  the  video  camera  at  Reynolds  numbers 
ranging  from  1 140  to  2240  with  no  imposed  unsteadiness.    Video  sequences  are  taken  at 
different  Reynolds  numbers  spaced  apart  by  approximately  50  units.  The  purpose  of 
this  flow  visualization  is  to  record  transition  events  with  no  imposed  unsteadiness. 
Still  photographs  are  shown  in  Figures  83-94. 
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IV.   RESULTS  AND  DISCUSSIONS 

The  coordinate  system  is  right  handed,    x  is  the  streamwise  distance  measured 
from  the  exit  end  of  the  nozzle,  y  is  the  vertical  distance  measured  from  the  bottom  wall 
of  the  channel,  and  z  is  the  spanwise  distance  measured  from  the  channel  centerline. 
The  positive  z  distances  are  to  the  right  of  the  channel  centerline  (looking  downstream). 

The  Reynolds  number.  Re,  is  based  on  the  bulk  flow  velocity,  vtxik.  and  the  channel 
height  d,  and  is  given  as: 

Re  =  vbulk  *  d/v 
Where  v  is  the  kinematic  viscosity.  The  Strouhal  number,  Str,  is  based  on  the  imposed 
unsteadiness  frequency,  fosc  .  bulk  flow  velocity,  and  channel  height.  The  Strouhal 
number  is  given  as: 

Str  =  2;:*f03c  *d/vbiik 
Reynolds  numbers  ranged  from  1 100  to  3800  in  the  present  study.  The  imposed 
oscillation  frequency  is  2.0  Hz.   Strouhal  numbers  varied  from  0.0351  to  0. 121. 

A.      CHANNEL  FLOW  QUALDTICATION 

For  fully  developed  laminar  channel  flow,  the  mean  velocity  is  given  by  as: 
U  /L'max^.Oly/d-fy/d)2] 
Where  u  is  the  local  mean  velocity  and  U  max  is  centerline  velocity  [Ref  ll:p.  280-281]. 
Figure  5  shows  the  velocity  profile  taken  at  x/d  =  300  at  a  Re  of  1 103  (vtnJk=  1.32  m/s) 
The  velocity  profile  (normalized  with  respect  to  U  rt]a^ )  is  parabolic  and  is  in  good 
agreement  with  the  above  equation.    For  this  situation,  the  maximum  centerline 
velocity  is  1.774  m/s  giving  a  ratio  of  v^ulk  /U  max  of  0.774. 
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B.  EFFECT  OF  2.0  HZ  IMPOSED  UNSTEADINESS  ON  LOCAL  MEAN 
VELOCnY  AT  y/d  =  0.88 

Figure  6  shows  variation  of  the  mean  velocity  at  y/d  =  0.88  with  and  without  the 
imposed  unsteadiness  for  varying  Re.  The  variation  of  the  bulk  mean  velocity  is  also 
shown  in  the  plot.   The  mean  velocity  shows  somewhat  non-linear  behaviour  for 
Reynolds  numbers  from  2000  to  2 100.  The  mean  velocity  distribution  is  not  changed  by 
imposed  bulk  flow  unsteadiness. 

C.  EFFECT  OF  2.0  HZ  IMPOSED  UNSTEADINESS  ON  RMS  VELOCITY  AT  y/d  =  0.88 


Figure  7  shows  the  variation  of  the  time-averaged  RMS  velocity  -\J  u'    normalized 
with  respect  to  the  time  averaged  velocity.  U  at  y/d  of  0.88  and  over  Reynolds  numbers 
from  1 100  to  3800.   Data  are  given  with  and  without  unsteadiness  at  2.0  Hz  at  twenty 
one  different  Reynolds  numbers. 

For  Reynolds  numbers  below  1400,  no  effect  of  the  imposed  unsteadiness  is  observed 
in  the  normalized  RMS  values.  With  the  2.0  Hz  unsteadiness,  significant  increases  in 
the  RMS  velocity  values  occur  for  Reynolds  numbers  from  1400  to  2005.  The  maximum 
increase  in  the  normalized  RMS  velocity  is  about  31.5%  over  that  without  the 
unsteadiness.  At  higher  Reynolds  numbers,  greater  than  2210,  however,  the  plot  shows 
that  normalized  RMS  values  decrease  with  2.0  Hz  unsteadiness.  This  reduction  is 
observed  for  Reynolds  numbers  up  to  2400.  Here,  the  maximum  decrease  in  normalized 
RMS  velocity  with  the  imposed  unsteadiness  is  about  5.4%  and  occurs  at  a  Reynolds 
number  of  2302.  For  Reynolds  numbers  greater  than  2400,  the  values  of  normalized 
RMS  are  nearly  the  same  with  and  without  the  imposed  unsteadiness. 
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D.      EFFECT  OF  REYNOLDS  NUMBER  ON  PHASE-AVERAGED  VELOCITY 

Figures  8  through  35  show  phase-averaged  velocity  traces  at  five  different  Reynolds 
numbers.  The  pase-averaged  quantities  are  plotted  against  the  phase  angle  of  the 
imposed  unsteadiness  for  Re  of  1 103,  1579.  2005,  2548.  and  3660. 

Figures  8  through  17  are  the  phase-averaged  quantities  with  no  imposed 
unsteadiness.   Figures  8,10,12,14.  and  16  show  phase-  averaged  velocity  0  (which 
equals  U   in  these  cases),  and  Figures  9,11,13,15,  and  17  show  the  plots  of  phase- 
averaged  RMS  velocity  normalized  with  respect  to  the  temporal  mean  U  .  For  these 
cases  phase  averaging  is  done  over  580  flow  oscillations.  The  variations  with  phase 
angle  in  Figures  8  through  1 7  are  quite  small  and  due  to  data  scatter. 

Figures  18  through  35  show  phase-averaged  velocity  traces  when  2.0  Hz 
unsteadiness  is  imposed  on  the  flow.  With  imposed  unsteadiness  at  this  frequency, 
statistically  stable  phase-averaged  results  are  obtained  after  580  cycles  for  Reynolds 
numbers  of  1200  and  below,  after  2320  cycles  for  Reynolds  numbers  up  to  1600,  and 
after  4640  cycles  for  higher  Reynolds  numbers  where  the  fluctuation  levels  are 
significant  with  respect  to  the  local  mean  velocity.   The  effect  of  multi-run  averaging  is 
demonstrated  in  Figures  20  through  35.   For  example.  Figures  24  and  25  show  phase- 
averaged  quantities  after  580  flow  oscillations,  and  Figures  26  and  27  show  the  same 
quantities  after  4640  flow  oscillations  for  Re  2005.  The  curves  in  figures  26  and  27  are 
smoother  indicating  improved  averages  of  0  and  normalized  RMS  values  at  particular 
phase  angles.  Figures  20  through  23  are  the  corresponding  plots  for  Re  1579,  Figures  28 
through  31  for  Re  2548,  and  Figures  32  through  35  for  Re  3660. 

Normalized  RMS  values  remain  relatively  constant  with  phase  angle  for  all  of  the 
cases  surveyed  with  the  exception  of  Re  of  2005.  Here,  as  shown  in  Figure  27,  an 
increase  in  the  normalized  RMS  velocity  is  seen  between  phase  angles  of  108  and  252 
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degrees.  The  maximum  value  Is  approximately  17%  greater  than  the  temporal  RMS 
velocity  and  occurs  at  a  phase  angle  of  about  200  degrees.  The  maximum  value  of  0   for 
this  case  is  seen  in  figure  26  at  about  235  degrees.  This  behavior  occurs  at  the  same 
Reynolds  number  where  the  time-averaged  RMS  velocity  shows  the  largest  increase  due 
to  imposed  bulk  unsteadiness. 

E.     MEAN  VELOCITY  PROFILES 

Figures  36  through  42  show  mean  velocity  profiles  taken  at  x/d  =  300.    y/d  ranges 
from  0.2  to  0.9  in  these  plots.  No  measurements  are  given  for  a  y/d  of  0. 1  because  of  the 
geometry  of  the  hot-wire  probe.   For  all  of  the  profiles,  mean  velocity  is  plotted  versus 
y/d  and  normalized  with  respect  to  the  maximum  mean  velocity. 

Figure  36  shows  the  profile  development  with  Reynolds  number  with  no  imposed 
unsteadiness.  Profiles  are  shown  for  Re  =  1103,  1579,  2005,  2548,  and  3660.  The  profile 
is  parabolic  at  low  Re  when  the  flow  is  fully  laminar.  At  the  highest  two  Re,  the  profiles 
show  fully  turbulent  behavior.   The  profile  is  distorted  for  the  transitional  Reynolds 
number  of  2005. 

Figure  37  shows  the  profile  development  with  Reynolds  number  with  imposed 
unsteadiness  at  2.0  Hz.  The  imposed  unsteadiness  has  no  observed  effect  on  the  mean 
velocity  profiles.  This  is  further  evidenced  in  Figures  38  through  42,  where  the  mean 
velocity  profiles  are  plotted,  both  with  and  without  imposed  unsteadiness,  for  each 
Reynolds  number. 
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F.  RMS  VELOCITY  PROFILES 

Figures  43  through  49  show  the  development  of  RMS  velocity  profiles  with 
Reynolds  number.   In  Figures  45  through  49,  RMS  velocities  are  normalized  with 
respect  to  maximum  observed  mean  velocities  and  plotted  against  y/d.  Again,  the 
symmetry  of  the  profiles  is  distorted  because  data  are  not  given  for  a  y/d  of  0. 1 . 

Figure  43  shows  dimensional  RMS  velocity  profiles  for  the  five  Reynolds  numbers 
when  no  bulk  unsteadiness  is  imposed  on  the  flow.   Corresponding  profiles  with 
imposed  unsteadiness  at  2.0  Hz.  are  shown  in  Figure  44.   In  Figures  45  through  49  plots 
both  with  and  without  imposed  unsteadiness  are  given  for  each  Reynolds  number.   For 
Re  =  2005  and  lower,  the  imposed  unsteadiness  results  in  significant  increases  in 
magnitudes  of  normalized  RMS  velocity  values.  At  Re  2005  values  range  from  0.04  to 
0.11  with  imposed  unsteadiness  compared  to  0.03  to  0.08  without,  and  the  relative 
shapes  of  the  two  profiles  are  nearly  the  same,  as  shown  in  Figure  47.  This  is  not  the 
case  for  Reynolds  numbers  1 103  and  1579. 

Figures  48  and  49  present  profiles  for  Reynolds  numbers  2548  and  3660.  At  these 
Reynolds  numbers,  the  imposed  unsteadiness  seems  to  have  little  effect  on  RMS 
velocity  distributions. 

G.  VARIATION  OF  PHASE-  AVERAGED  QUANTITIES  WITH  REYNOLDS  NUMBER 
AND  y/d 

Figures  50  through  82  show  the  variations  of  phase-averaged  quantities  u  and  u  as 
y/d  and  Reynolds  number  are  changed. 

Figures  50  through  8 1  show  0  ,  and  0  '  (normalized  with  respect  to  the  temporal 
mean)  plotted  versus  phase  angle  at  Re  2005  for  y/d  from  0.2  to  0.9.  The  figures  are 
arranged  to  show  data  as  y/d  decreases. 
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Figures  50  through  57  are  the  plots  of  0  versus  phase  angle  for  flow  with  no 
imposed  unsteadiness.  The  corresponding  0  '  plots  are  shown  in  Figures  58  through  65. 
Phase-averaged  RMS  U  '  decreases  from  y/d  0.9  to  0.7,  and  then  becomes  relatively 
constant  at  y/d  of  0.6,  0.5,  and  0.4.  For  flow  with  no  imposed  unsteadiness,  the  0  ' 
distributions  are  relatively  constant  with  phase  angle  for  each  y/d  location. 

Figures  66  through  73  show  0  distributions  versus  phase  angle  for  Re  of  2005  with 
imposed  unsteadiness  at  2.0  Hz  and  at  y/d  locations  0.9  to  0.2.  Maximum  values  of  0 
are  observed  at  phase  angles  between  about  210  degrees  and  240  degrees.   Minimum 
values  occur  at  phase  angles  between  approximately  30  and  75  degrees.   Near  the  wall 
the  curves  show  considerable  fluctuations  as  shown  in  Figures  66,  67,  and  73.  Figures 
68  through  72  show  smoother  curves  as  the  channel  centerline  is  approached.  Most  u 
distributions  are  nearly  sinusoidal  in  shape,  as  shown  in  Figure  70. 

Figures  74  through  81  show  plots  of  normalized  U  '  versus  phase  angle  with 
imposed  unsteadiness  at  2.0  Hz.  The  previously  described  increase  in  phase-averaged 
RMS  velocity  at  this  transitional  Reynolds  number  is  again  evident  between  phase 
angles  of  108  and  252  degrees.  This  increase  is  less  dramatic  at  y/d  locations  near  the 
channel  centerline.    In  fact,  the  distribution  is  relatively  constant  with  phase  angle  at 
y/d  of  0.5,  as  shown  in  Figure  78.  The  most  noticeable  increases  are  observed  at  y/d 
locations  of  0.9,  0.8,  and  0.2,  as  shown  in  Figures  74,  75,  and  81,  respectively. 

In  Figure  82,  peak  to  peak  phase-averaged  velocities  0  (normalized  with  respect  to 
the  local  mean  velocity)  are  shown.   Data  for  2.0  Hz  imposed  unsteadiness  are  given  for 
Re  of  1 103,  1579,  2005,  2548,  and  3660.  The  y/d  locations  vary  from  0.5  to  0.9.  For  the 
laminar  cases.  Re  of  1 103  and  1579,  the  peak  to  peak  velocity  shows  a  steady  increase 
from  6%  at  the  channel  centerline  to  about  12%  near  the  top  wall  (y/d=0.9).  For  the 
intermediate  Re  of  2005,  the  variation  is  greater,  11%  at  the  centerline  and  24%  near 
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the  wall.  For  Re  2548.  the  variation  at  the  centerline  of  the  channel  is  the  same  as  that 
for  Re  2005.   However,  at  the  remaining  y/d  locations  the  respective  values  are 
significantly  less  than  those  for  Re  2005.  At  y/d  of  0.9  the  variation  for  Re  2548  is  7% 
less  than  that  for  Re  2005  at  this  location.  The  shape  of  the  curve  for  Re  2548  is  similar 
to  the  curves  observed  for  Reynolds  numbers  1 103  and  1579.  For  Re  3660  the 
distribution  is  nearly  constant  at  about  12%  from  y/d  0.5  to  0.8.  Then,  peak  to  peak 
magnitude  increases  sharply,  reaching  a  maximum  value  of  25%  at  y/d  of  0.9. 

H.      FLOW  VISUALIZATION 

Figures  83  through  94  show  flow  visualization  results  with  no  imposed 
unsteadiness  at  Reynolds  numbers  ranging  from  1 136  to  2239.   Flow  is  from  the  right  to 
the  left  of  the  photographs  in  Figures  83  through  92.  In  Figures  93  and  94,  flow  is  from 
the  bottom  to  the  top  for  each  photograph. 

Figures  83  through  92  show  flow  structures  occuring  during  the  transition  process 
from  laminar  to  turbulent  flow.   Laminar  flow  at  Re  1 136  is  shown  in  Figure  83.   Figure 
84  shows  the  presence  of  Tollmien-Schlichting  waves  at  the  higher  Reynolds  number. 
Re  1328.  The  emergence  of  vortical  ribbon-like  structures  (hereafter  referred  to  as 
ribbons)  at  the  left  side  of  the  channel  is  seen  at  Re  1377,  in  Figure  85.  As  Re  increases, 
these  ribbons  are  observed  to  occupy  a  greater  percentage  of  the  width  of  the  channel, 
from  left  to  right.  This  is  shown  in  Figures  86  and  87  for  Re  of  1395  and  1496, 
respectively.  Also  noticeable  at  Re  1496  is  the  presence  of  turbulent  spots  which  appear 
as  dark  areas  in  Figure  87.  At  higher  Reynolds  numbers,  turbulent  spots  are  larger  and 
more  abundant.  This  is  shown  in  Figures  88  through  91,  which  correspond  to  Re  of 
1529,  1713,  and  1854,  respectively.  Figure  92  shows  a  predominantly  turbulent  flow  at 
Re  2239.  Turbulent  ribbons  which  are  hazy  and  ill-defined  compared  to  those  observed 
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at  lower  Reynolds  numbers  are  characteristic  at  this  Reynolds  number,  as  shown 
in  the  figure. 

Figures  93  and  94  show  a  time  sequence  which  depicts  flow  behavior  at  a  sub- 
critical  Reynolds  number  of  1412.  The  time  interval  between  each  of  the  twelve  frames 
of  the  sequence  (ordered  from  top  to  bottom  of  the  page,  left  to  right)  is  approximately 
one  fifteenth  of  a  second.  Therefore,  the  total  time  sequence  represents  flow  history  of 
approximately  800  milliseconds.  The  sequence  clearly  illustrates  the  order  and 
duration  of  transition  phenomena  with  no  imposed  unsteadiness,  as  soon  as  the  smoke 
is  released  from  the  wire.  The  pictures  suggest  that  the  ribbons  are  considerably 
elongated  in  the  streamwise  direction  and  occupy  half  of  the  channel  width.   Some  of 
these  ribbons  terminate  as  dark  patches. 
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V.  CONCLUSIONS 

Imposed  unsteadiness  at  2.0  Hz  has  no  discernible  effect  on  mean  velocity  profiles  for 
Re  from  1 103  to  3660.   However,  in  regard  to  time-averaged  fluctuating  velocity 
magnitudes,  the  imposed  unsteadiness  exhibits  both  a  destabilizing  effect  and  a 
stabilizing  effect  in  the  range  of  Reynolds  numbers  and  Strouhal  numbers  studied.  A 
destabilizing  of  the  flow  is  seen  for  Re  from  1400  to  2200,  and  Str  from  0.061  to  0.095, 
whereas  a  stabilizing  of  the  flow  is  seen  for  Re  from  2200  to  2400  and  Str  from  0.056  to 
0.061. 

The  imposed  unsteadiness  affects  both  time-averaged  and  phase-averaged  RMS 
values.  This  is  especially  evident  for  the  transition  Reynolds  number  of  2005.   Here, 
when  the  unsteadiness  is  imposed,  the  phase-averaged  RMS  velocity  is  a  maximum  at  a 
phase  angle  slightly  ahead  of  the  phase  angle  where  the  maximum  phase-averaged 
velocity  0   is  observed.   The  most  significant  variations  of  the  phase-averaged  velocity 
across  the  span  of  the  channel  also  occurs  at  this  Reynolds  number. 

The  various  transition  flow  structures,  which  are  predicted  by  numerical  simulations 
and  have  been  observed  experimentally  during  previous  flow  visualization 
investigations,  are  confirmed  in  this  study.    Tollmien-Schlichting  waves,  followed  by 
vortical  ribbon-like  structures,  and  then,  turbulent  spots  are  the  dominant  phenomena 
as  Reynolds  number  is  increased  with  no  imposed  unsteadiness.  Time  sequence 
photography  at  a  sub-critical  Reynolds  number  provides  insight  into  the  emergence  of 
the  ribbon-like  vortical  structure  in  the  channel,  with  subsequent  growth  in  both  the 
spanwise  and  streamwise  directions.   Dark  patches  are  sometimes  observed  following 
the  ribbon  structures,  suggesting  the  onset  of  short  duration  turbulent  spots. 
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While  the  above  results  proved  repeatable  during  this  study  more  investigation  is 
warranted,  especially  in  the  range  of  Reynolds  numbers  where  the  imposed 
unsteadiness  seems  to  stabilize  the  flow,  i.e..  Re  2200  to  2400. 


27 


APPENDIX  A 


SOFTWARE     DIRECTORY 


1 .  HWCAL  :  This  program  determines  the  constants  for  the  King's  Law  calibration  of 
the  hot-wire.   The  program   also  provides  a  polynomial  fit  of  the  calibration  data. 

2.  HOTWTREPAV  :  This  program  is  used  to  read  the  data  stored  in  the  A/D  buffer  of  the 
high  speed  data  acquisition  system  and  stores  the  information  on  micro  diskettes. 
Manual  inputs  are:  triggering  frequency,  hot-wire  DC  voltage  (ungained),  oscillation 
frequency  (flow  blockage),  bulk  velocity  and  date  and  time  of  run. 

3.  MAX-MIN  :  This  program  calculates  instantaneous  and  phase  averaged  velocities. 
Initially  a  look  up  table  is  created.    Here,  effective  velocities  are  calculated  from  the 
effective  voltage  values  and  stored  for  follow  on  calculations.   The  hot-wire  calibration 
constants  obtained  from  HWCAL,  and  the  amplifier  gain  are  incorporated  into  these 
calculations.    The  velocity  calibration  is  given  by  the  equation: 

uerr=k(Een-2-Eo2)1/N 

where  k  is  the  proportionality  constant,  Eeflis  the  effective  voltage,  and  E0  is  the 
reference  voltage  at  no  flow.   N  is  a  constant  value  of  0.45  for  moderate  Reynolds 
numbers.  Once  the  look  up  table  is  created  the  program  reads  the  instantaneous  voltage 
values  from  the  data  file  and  converts  them  to  instantaneous  velocities.  At  this  point  a 
plot  of  the  instantaneous  velocity  versus  time  can  be  generated  for  any  of  580  cycles 
(The  first  60  cycles  are  discarded  to  allow  for  flow  stabilization  as  the  unsteady  device  is 
started.) 
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Next  the  program  phase  averages  the  580  cycles,  and  velocity  versus  phase  angle  (of 
the  flow  blockage)  plots  are  available.  Two  plots  are  available  from  the  averaged 
values;  u  versus  phase  angle  and  V  u1     divided  by  U,  versus  phase  angle,  where  u  is  the 
phase  averaged  velocity,    Vu      is  the  phase  averaged  root  mean  squared  velocity,  and 
U  is  the  average  velocity.  In  the  case  where  there  is  no  imposed  unsteadiness,  the  phase 
averaged  velocity,  u.  is  equal  to  the  time  averaged  velocity  ,  U   since  there  is  no 
superimposed  phase  average  velocity  component,  Q  • 

The  program  can  also  be  used  to  process  a  set  of  data  when  multi-run  averaging  is 
desired.  In  this  case  the  averaged  values  are  written  to  and  stored  in  the  computer  hard 
disk  memory.  Once  each  of  the  desired  data  sets  has  been  averaged  and  stored  in  the 
computer  memory,  they  are  copied  to  micro  diskettes. 

4.  AVGPLOT  :  The  program  AVGPLOT  is  used  to  average  any  number  of  data  sets.  A 
maximum  of  eight  sets  of  data  (or  4,640  cycles)  were  averaged  for  this  study.  The 
program  AVGPLOT  provides  the  same  graphs  as  those  available  from  MAX_MIN.   Phase 
averaged  data  sets  are  copied  to  micro  diskettes  prior  to  purging  them  from  the 
computers  memory. 
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APPENDIX  B 


FIGURES 
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Figure  1.   Straight  Channel 
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Figure  2.   Unsteady  Device 
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SCHEMATIC  OF  FLOW  MEASUREMENT  EQUIPMENT 


Figure   3 
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Figure   83.      Re  =   1136 
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Figure   84.      Re   =    1328 
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Figure   85.      Re  =   1377 
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Figure    86.      Re   =    1395 
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Figure   87.      Re   =    1496 
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Figure   88.      Re  =   1529 
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Figure   89.      Re  =   1713 
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Figure    90.      Re   =    1713 
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Figure  91.   Re  =  1854 
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Figure    92.      Re   =   2239 
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Figure  93.   Time  Sequence  at  Re  =  1412 
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Figure  94.   Time  Sequence  at  Re  =  1412,  continued. 
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